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Viscous  hypersonic  f  low  at  high  Reynolds  numbers  R  (R  =  Pouo^c/r‘) 

over  a  blunt-nosed  body  is  considered  on  the  basis  of  Navier-Stokes 
equations.  An  asymptotic  solution  of  the  Navier-Stokes  equations  Is 
obtained  for  the  case  where  substantial  viscous  effects  are  present  in 
the  entire  Inner  layer. of  rthickness  (e  is  the  shock  density 

ratio).  This  means  that  the  thickness  of  the  boundary  and  the  Inner 
inviscid  layers  are  of  the  samooorder  of  magnitude'.  By  substituting 
the  nondimensional  coordinates  and  functions  of  the  viscous  flow 
region  in  the  Navier-Stokes  equations  and  taking  the  similarity 
/  5/2 

parameter  A  =  1/R^e  '  constant.  Into  account.,  a  system  of  equations 

is  derived.  A  distinctive  feature  of  the  boundary  value  problem 
considered  here  is  that  It.  contains  a  new  similarity  parameter  a .  It. 
is  demonstrated  that  the  asymptotic  solution  whore  A  tends  to  aero 
coincides  with  that  of  Fay  and  Riddel.  The  numerical  results  of 
computations  on  an  M-20  computer  for  various  at  three  different. 
temperature1  parameters  gu>(gu>  =  1  corresponds  to  adiabatic  wall), 
presented  in  tabular  form,  show  that.  at.  A  =  1  the  heat  flux  Is  17-23$ 
higher  than  In  the  case  of  an  ordinary  boundary  layer.  The  value 
A  ~  1  corresponds,  for’  example,  to  flight  conditions  at  M-20, 

+4 

Re  =  10  ,  and  e  =  0.1.  Orig.  art.  has:  13  formulas. 
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CALCULATION  OF  HEAT  TRANSFER  ON  FRONTAL  SURFACE 
OF  BLUFF  BODY  IN  HYPERSONIC  FLOW 

Yu.N.  Yermak,  V.Ya.,  Neyland 

(Moscow) 


The  present  paper  considers  hypersonic  flow  of  a  viscous  gas 
past  the  nose  of  a  bluff  body  at  high  Reynolds  numbers  /?(/?  = 

=  poueUo/r,  where  the  subscript  0  is  attached  to  the  values  of 
quantities  in  the  oncoming  flow) .  A  rather  large  number  of  works, 
among  which  we  may  cite  [1-3],  has  been  devoted  to  this  problem. 
The  purpose  of  the  investigations  usually  consists  in  determining 
heat  flow  at  the  critical  point  of  the  body.  These  studies  have 
been  made  within  the  framework  of  boundary-layer  theory  in  the 
first  [1,  2]  or  second  [3]  approximation.  When  we  remember  that 
at  hypersonic  flow  velocity,  the  ratio  of  the  gas  densities  in 
front  of  and  behind  the  compression  shock  e  *=  p«/pi  is  small,  the 
solution  of  the  problem  can  be  approached  by  the  so-called  New¬ 
tonian  theory,  within  whose  framework  we  seek  an  asymptotic  solu¬ 
tion  to  the  Navier-Stokes  equations  (see,  for  example,  [3])  as 
Wo  e  -►  0  and  /?-*•».  In  first  approximation  in  this  theory, 

an  inviscid  shock  layer  with  a  thickness  of  the  order  of  er  is 
formed  between  the  body  and  the  shock  wave.  The  flow  in  it  is  de¬ 
scribed  around  the  critical  point  by  the  Euler  equations,  in  which 
the  term  taking  account  of  the  longitudinal  pressure  gradient  van¬ 
ishes.  To  obtain  these  equations.  It  is  necessary  to  introduce  di¬ 
mensionless  coordinates  with  consideration  of  shock-layer  thick¬ 
ness  and  the  longitudinal  dimension  of  the  minimum  Influence  re¬ 
gion  near  the  compression  shock: 


X  =  rc'l >x,  y  «=  rey  (1) 

where  r  is  the  radius  of  curvature  of  the  body’s  contour  and  x,  y 
are  orthogonal  curvilinear  coordinates  bound  to  the  surface  of  the 
body  in  the  direction  of  the  normal  to  it.  The  expansions  for  the 
flow  function  in  the  neighborhood  of  the  critical  point  can  be 
presented  in  the  form 

u  «=>  uoe'/i (u(f)5  +  i>  =  u,e[r(jr)  +  ...) 

(2) 

p  ~  />  *=  i>ou»: !/'(;/>  +  •••].  /» «=  +  .••) 

Then  we  obtain  in  first  approximation  from  the  Navior-Stokes 
equations 
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2a+  v'  =  0,  u:  +  va'  =  0  (  3  ) 

Here  we  have  dropped  the  prime  on  the  dimensionless  quanti¬ 
ties.  Considering  the  boundary  conditions  at  the  compression 
shock,  we  can  obtain 

u  =  y~  +  (ij) 

since  on  the  wave  u  =  -v  *  1  for  y  =  y 2.  Near  the  body's  surface 
v  -*•  0,  u  ■*  0. 

The  inertial  terms  in  the  momentum  equation  in  a  layer  with 
a  thickness  of  the  order  of  re*'2  acquire  the  same  order  of  magni¬ 
tude  as  the  dropped  term  with  the  longitudinal  pressure  gradient. 
Hence  to  obtain  the  uniform  exact  solution  in  this  region  as  well, 
it  is  necessary  to  introduce  new  dimensionless  coordinates  and 
flow  functions: 

*  *“  rt'ltx.  jr  bi  a  «=  u,t[U(Y)x  +  ...1  o  «=  o«e1[V(y)  +  ...]  (5) 

Such  a  layer  was  examined  in  [*J]  in  solving  the  problem  of 
hypersonic  flow  of  an  inviscid  gas  past  a  flat  plate  set  trans¬ 
verse  to  the  oncoming  flow.  After  substitution  of  (*3)  in  the 
Navier-Stokes  equations,  we  obtain,  in  first  approximation. 


2  U  +  V  —  0. 

U*+VU'=>  2 

(6) 

Integrating  (6),  we  obtain 

=  iY+yfc)ic, 

V  «=  C(2  —  IP) 

(7) 

Use  of  the  familiar  principle  of  splicing  the  external  and 
internal  asymptotic  expansions  enables  us  to  determine  the  con¬ 
stants:  C  =  1  and  tjz  =  1  +  0{e'/i).  Thus  we  obtain  the  velocity  gradient 
at  the  body's  surface  as  u^/Tz/v.  This  value  is  also  obtained  in 
the  conventional  Newtonian  theory.  It  is  used  as  the  boundary  con¬ 
dition  for  the  boundary  layer,  whose  equations  are  obtained  by 
asymptotic  expansion  of  the  Navier-Stokes  equations  in  a  region 
with  a  thickness  of  order  R~%  immediately  adjacent  to  the  body's 
surface. 

Thus,  In  the  conventional  scheme  for  finding  the  asymptotic 
solution  of  the  Navier-Stokes  equatl  on:’. ,  the  limit  transition  n  -► 

-*•  0,  Wo  -*■  °°,  R  -*•  Is  accomplished  in  such  a  way  that  the  entire 
shock  layer  will  decompose  Into  two  inviscid  regions  and  a  bound¬ 
ary  layer.  This  limit  transition  describes  flow  conditions  well 
for  Mach  number  Mo  and  r  moderate  and  R  large. 

The  present  work  examines  another  type  of  limit  transition, 

-*•  00 ,  e  0 5  R  -*•  °°,  in  which  viscous  effects  remain  essential 
throughout  an  inner  layer  with  thickness  e3/2r.  This  means  that 
the  thickness  of  the  boundary  layer  and  the  internal  inviscid 
layer  are  of  the  same  order  of  magnitude. 

Let  us  introduce  the  parameter  ih  =  /tin / no.  In  the  limit  transi¬ 
tion  under  examination,  the  parameter  a  =  1  remains  finite  as 

Mo  -►  °°,  e  -*•  0  and  R\  -*•  °°.  It  characterizes  the  ratio  of  the  thick- 
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nesses  of  the  boundary  and  inner  inviscid  layers.  As  estimates 
show,  the  thickness  of  the  compression  shock  remains  negligibly 
small,  i.e.,  passage  across  the  compression  shock  is  described, 
as  before,  by  the  Hugoniot  solution.  The  asymptotic  expansions 
and  the  solutions  retain  their  previous  form  ( 1 ) — ( 3 )  in  a  layer 
whose  thickness  is  of  the  order  of  er. 

The  dimensionless  coordinates  and  functions  in  the  viscous 
flow  region  have  the  form  (5).  Substituting  them  into  the  Navier- 
Stokes  equations  and  considering  the  condition  —  const,  we  obtain, 
in  first  approximation,  the  following  equations,  which  have  been 
expanded  in  x: 

2pU+  (pV)'~0,  pVf'  «■  A[(|i/ a) cl',  p(P*+  VIT)  -  2A(ntT)'  (8) 

The  boundary  conditions  for  System  (8)  are  known  only  at  the 
wall,  V  =  v  =  D,  and  the  missing  ones  are  obtained  by  splicing 
with  Solution  ( t) .  Introducing  the  stream  function  and  Dorodni- 
tsyn's  variables,  we  transform  the  system  to 

Y 

A  W)'-r2//'-!-2/P-frF«0.  A  fhV/«)  (.'T  l" 0  (W  =  .  r|  --=i  \  pdY j  (9) 


Here  differentiation  is  with  respect  t'o  the  Dorodnitsyn  vari¬ 
able.  The  boundary  conditions  at  the  wall  have  the  usual  form 

/( 0)“/'(0)*,0;  p(0)=ff„  (10) 

The  missing  boundary  condition  for  the  momentum  equation, 
which  is  obtained  by  splicing  the  solutions  of  System  (9)  with 
(;4),  takes  the  form 

~  «=*  ( Ui:  /  0;i)  „  .0  1,  or  /"(«.)  =  1  ( 11 ) 

The  external  boundary  condition  for  the  energy  equation  has 
remained  as  before: 

*.’(<«)  «  1  —  (12) 

The  boundary -value  problem  (9)-(12)  differs  from  that  con¬ 
sidered  earlier  in  [1,  2]  and  others  in  that  the  new  similarity 
parameter  A  has  appeared;  its  physical  sense  Is  indicated  above. 
The  external  boundary  condition  for  the  momentum  equation  has  also 
changed.  It  is  evident  from  physical  considerations  that  the 
asymptotic  solution  (9) -(12)  must  be  identical  to  the  solution  ob¬ 
tained  In  [1,  2]  and  other  papers  as  A  -*•  0 .  It  is  not  difficult 
to  show  that  this  is  Indeed  the  case,  since  the  asymptotic  solu¬ 
tion  (9)— (12)  takes  the  following  form  as  A  •+  0: 

/"-»<?" /Va7  tf-^h’iyiC  (13) 

where  and  h  are  solutions  to  the  boundary-value  problem 
(.\y ') '  +  +  2  / ,»  -  («F')»  -  0,  [  (/v  /  cr)  +  2<p/.'  0 

qi(0)  Cf'(0)  0,  //  (0)  hw  H ,  fp'(oo)  m  y*2,  A(oo)  1 
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Equations  (9)  were  solved  on  the  M-20  electronic  computer 
[EC]  (3BH)  with  boundary  conditions  (10),  (11)  and  (12).  We  pre¬ 
sent  the  results  of  calculations  for  various  A  and  three  values 
of  the  temperature  factor  gu  (the  value  gw  =  1  corresponds  to  an 

adiabatic  wall) . 
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W* 

© 

H 

•0.88, 

4-0.41 

A -10“*  0.37 

10*0.11 

10*0.22 

0.11 

0.43 

0.51 

C= 14.82 

8.76 

6.31 

3.11 

1.80 

1.70 

14.3# 

8.39 

5.40 

2.67 

1.33 

1.23 

3.92 

2.79 

1.31 

0.69 

0.64 

3.90 

2.765 

1.24 

0.62 

0.57 

A  -  0.63 

0.79 

0.00  1.00 

3.86 

11.33 

4-1.58 

1.46 

1.41  1.36 

0.96 

o.so 

%//S-1.07 

0.98 

0.92  0.88 

0.447 

0.262 

4-0.58 

0.53 

0.50  0.48 

0.27 

0.18 

-0.517 

0.45 

0.43  0.41 

0.209 

0.13 

?-=io-*. 

4  =0.24, 

*^=0.13 

A -10*0.16  10*0.24 

10*0.17  10"*0.19  10*0.51 

0.10 

4-63.36 

5.25 

2.09 

1.98 

1.26 

0.92 

4/V7S-61.53 

4.85 

1.80 

1.72 

1.05 

0.76 

4-34.51 

2.80 

1.09 

1.02 

0.63 

0.45 

HJV*-  33.33 

2.61 

1.01  0.93 

0.57 

0.41 

A  =  0.19 

0.48 

0.R8  1.00 

2.39 

4.48 

4-0-72 

0.51 

0.42  0.40 

0.32 

0.28 

~4/VA«0.55 

0.347 

0.256  0.24 

0.15 

0.11 

4-0.34 

0.22 

0.17  0.18 

0.11 

1(7*0.89 

-4/VS-0.298 

0.187 

0.138  0.18 

0.84.1(7* 

0.01.1(7* 

k'm  -1. 

*w=*o,  • 

4=2.200 

A-0.8M0* 

0.44-10*  0.20-10* 

0.40-10*  0.13  0.25 

4 « 77.91 

33.79  10.24 

10.91 

0.82  5.10 

~4//A=  77.40 

33.2  15.63 

10.25 

0.13  4.40 

A  =0.59  0.7 

0.814  1.00 

1.9' 3 

10 

4-3-5' 

p  3. 346 

3.153  2.92 

2.30 

1.499 

• 

~4//5-2.87  2. 63 

2.44  2.200 

1.50 

0.09 

It  was  assumed  in  all  calculations  that  o  =  0.7^  and  that  the 
variation  of  viscosity  p  'v  gu,  where  u>  =  0.76. 

It  can  be  seen  that  for  A  =  1,  the  heat  flow  is  17-23?  larger 
than  the  corresponding  value  for  the  ordinary  boundary  layer.  The 
quantity  A  'v  1,  for  example,  corresponds  to  flight  regimen  with 
M  >  20,  e  *  0.1  and  Reynolds  numbers  R  of  the  order  of  10+‘’. 
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